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Abstract 

We present the state-of-the-art of the effective field theory computation of the MSSM 
Higgs mass, improving the existing ones by including extra threshold corrections. We show 
that, with this approach, the theoretical uncertainty is within 1 GeV in most of the rele¬ 
vant parameter space. We confirm the smaller value of the Higgs mass found in the EFT 
computations, which implies a slightly heavier SUSY scale. We study the large tan j3 region, 
finding that sbottom thresholds might relax the upper bound on the scale of SUSY. We 
present SusyHD, a fast computer code that computes the Higgs mass and its uncertainty 
for any SUSY scale, from the TeV to the Planck scale, even in Split SUSY, both in the DR 
and in the on-shell schemes. Finally, we apply our results to derive bounds on some well 
motivated SUSY models, in particular we show how the value of the Higgs mass allows to 
determine the complete spectrum in minimal gauge mediation. 
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1 Introduction 

With the recent discovery of the Higgs boson the last missing piece of the Standard Model (SM) has 
been unveiled and all the parameters of the theory measured. The success of the Standard Model 
in describing all the measured observables at colliders contrasts with the failure of the theory to 
explain some non-collider observations, such as Dark Matter, matter-antimatter asymmetry, etc. 

Among the various completions of the Standard Model proposed so far, supersymmetric 
(SUSY) theories remain the most attractive option. Not only they screen the electroweak (EW) 
scale from ultraviolet (UV) sensitivity to new physics thresholds but they successfully predict the 
unihcation of gauge couplings and may provide with a natural WIMP dark matter candidate. 

On the other hand, indirect hints for a small hierarchy between the scales of electroweak and 
supersymmetry restoration (e.g. flavor observables, searches for EDMs, LEP bounds on the Higgs 
mass), have found stronger support from the recent discovery of a moderately heavy SM-like 
Higgs [1,2] and from the absence of any evidence of superpartners in the hrst LHC runs at 7 and 
8 TeV (see e.g. [3,4]). 
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While the scale of supersymmetry may still be low (there are still various arguments in favor of 
this scenario), hopefully within the reach of the next LHC run, it is fair to say that our conhdence 
in predicting the new physics scale based on naturalness arguments weakened substantially [5,6]. 
The suspicion that other mechanisms may explain the strength of the weak interactions is becoming 
stronger and alternative scenarios to low energy SUSY already exist [7,8]. It is thus useful to look 
for different (more-experimentally-driven) methods to infer the scale of the superpartners. One 
natural candidate is the value of the Higgs mass, which in supersymmetry is calculable in terms 
of the couplings and the soft SUSY breaking parameters. 

In particular, in the Minimal Supersymmetric Standard Model (MSSM) the tree-level Higgs 
mass is predicted to lie below the Z-boson mass up to quantum corrections logarithmically sensitive 
to the SUSY breaking scale. Therefore the measured value of the Higgs mass gives non-trivial 
constraints on the spectrum and couplings of the MSSM, allowing to shrink the allowed energy 
range for the superpartners. 

Given the logarithmic (in)sensitivity of the Higgs mass value to the SUSY scale, high precision 
is required in such calculation to reliably determine the allowed parameter space of the theory. Be¬ 
sides, the experimental value is now known with per mille accuracy ruh = 125.09(24) GeV [9]. The 
effort in the Higgs mass calculation has been remarkable, reaching the two-loop and in some cases 
the three-loop level, with different techniques and schemes, see e.g. [10-26]. Some of the compu¬ 
tations, however, are only valid for small SUSY breaking scales, where log-corrections do not need 
resummations; in fact currently available computer codes have a very limited range of applicabil¬ 
ity compared to the allowed parameter space. Moreover, different computations and/or computer 
codes disagree among themselves, in some cases substantially more than the expected/claimed 
level of uncertainties. 

Given the important role played by the Higgs mass in constraining supersymmetric models, 
the limitations of the existing codes and the disagreements in the literature, we felt the need to 
revisit the computation. We put special emphasis on the relevant parameter space to reproduce 
the experimental value of the Higgs mass, on the study of the uncertainties and on the possible 
origin of the differences with other methods. In this paper we recompute the Higgs mass in the 
MSSM using the effective held theory (EFT) approach, which allows to systematically resum large 
logarithms and to have arbitrary big hierarchies in the spectrum, exploiting the mass gap hinted 
by the largish value of the Higgs mass and the absence of new physics at the LHG. 

Our computation follows very closely the ones in [7,8,10,23-26], providing independent checks 
of such computations. We improved them in various ways. We added all the dominant SUSY 
threshold corrections including the contributions from bottom and tau sectors, which become 
important at large tan /3. In this way we provide the state of the art in the EFT calculation of 
the Higgs mass. We performed the computation in both the DR and the on-shell (OS) schemes, 
the latter has the advantages of ensuring the correct decoupling limits and keeping the theoretical 
errors under control in the whole parameter space. We point out that large logarithms arising 
from splitting the fermions from the scalar superpartners in split SUSY scenarios do not need 
resummation in the whole region of parameters space relevant for the observed Higgs mass. We 
also hnd that nih = 125 GeV may not necessarily bound the SUSY scale to he below 10^° GeV 
(and much below at large tan/?), but it might extend at arbitrarily high scales. Another outcome 
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of our computation is that, even at maximal stop mixing, the average stop mass is required to be 
above the TeV scale in order to reproduce the correct Higgs mass in the MSSM. We also performed 
a study of the various possible uncertainties, showing that for most of the parameter space we are 
dominated by the experimental ones. We identihed some of the sources of disagreement between 
existing computations/codes. 

We implemented the computation in a new computer code, SusyHD, which we make public 
[27] and which allows to reliably compute the MSSM Higgs mass (and its uncertainties) even when 
big hierarchies are present in the spectrum. Avoiding slow numerical integrations, the code is fast 
enough to be used to set the experimental value of the Higgs mass as a constraint on any other 
SUSY parameter. 

Finally we also explore the implications of the Higgs mass on two of the simplest SUSY 
scenarios: minimal gauge mediation (MGM) and anomaly mediation. In particular in the first 
case we show how the value of the Higgs mass allows to determine the complete spectrum of the 
superpartners. 

The paper is organized as follows. In section 2 we describe our computation of the Higgs mass 
in the effective held theory approach, we study the theoretical uncertainties and we compare our 
results with the existing ones. In section 3 we present the implications of our computation for the 
SUSY spectrum in different regimes, in particular we show the constraints from the Higgs mass in 
the parameter range relevant for SUSY searches at hadron colliders, we explored the region of very 
large tan/3 and we comment on the (non) importance of extra log resummation when the SUSY 
spectrum is split. In section 4 we briehy introduce SusyHD, a new code to compute the Higgs 
mass using the EFT technique. In section 5 we apply our results to two of the simplest SUSY 
models: minimal gauge mediation and anomaly mediation. We summarize the most signihcant 
results in the conclusions in sec. 6. Finally, in appendix A, we provide the explicit expressions 
for some of the SUSY thresholds computed in this work and more details about the conventions 
used in the text. The reader not interested in the technical details of the computation can look 
directly at secs. 3 and 5. 


2 The computation 

2.1 The Effective Field Theory technique 

Whenever a theory presents a gap in its energy spectrum effective held theory techniques become 
a very powerful tool. They exploit the hierarchy of scales to allow a perturbative expansion in 
powers of the energy gap. This simplihes the theory getting rid of irrelevant degrees of freedom 
and couplings. 

Applied to supersymmetry, as the scale of the superpartners is raised, the Standard Model 
becomes a better and better EFT, with corrections from higher dimensional operators decoupling 
fast, as powers of u/msusY, the ratio between the EW and the SUSY scale. At leading order in this 
expansion the presence of supersymmetry at low energy reduces to a boundary condition for the 
SM couplings evolved at the SUSY scale, where they have to match with the full supersymmetric 


3 



theory. 

From the bottom up the technique reduces to taking the measured SM couplings at low energy, 
evolving them up to the superpartner scale and matching them to the full supersymmetric theory 
living at high scales. The non-trivial relations between the couplings in the supersymmetric 
theory (in particular between the Higgs quartic, the gauge-Yukawa couplings and the soft terms) 
translates into a non-trivial condition on the soft SUSY parameters. Eqnivalently one can leave 
the physical Higgs mass as a free parameter to be determined as a fnnction of the UV SUSY 
parameters. Imposing the physical valne for the Higgs mass then gives the constraint. 

The nse of this technique in the compntation of the Higgs mass in the MSSM is qnite old [10], 
however its ntility in natural SUSY spectra was limited since corrections from higher dimensional 
operators could not be neglected in that case. These techniques became more popular with the 
advent of Split SUSY scenarios [7,8] and the recent LHC resnlts [24-26,28,29]. 

In the rest of the paper, nnless specihed otherwise, the gange couplings ( 71 , 2 , 3 , the Yukawa 
couplings yt,h,T and the Higgs qnartic conpling A are assnmed to be the SM ones in the MS scheme 
while the soft parameters (masses and trilinear couplings) are in the DR or OS schemes. In 
particular when we refer to onr DR or OS results it means that the soft masses are DR or OS 
while the couplings are always taken to be the SM ones in the MS scheme. 

Onr compntation is organized as follows: 

• The SM couplings (gauge, Yukawa and qnartic) in the MS scheme are extracted from the 
corresponding physical qnantities at the EW scale at fnll two-loop level [30]. In particular 
the matching between the top mass^ and the top Yukawa coupling is done using fnll two-loop 
thresholds plus the leading three-loop QCD one^ from [32]. 

• The conplings are then evolved from the weak scale to the snperpartner scale using the full 
three-loop renormalization gronp eqnations (RGE) for these conplings^ [30]. 

• At the SUSY scale the SM conplings are matched to those of the SUSY theory (converted 

from either DR or OS to the MS scheme) using the fnll one loop thresholds (from [26] 
and the corrections from appendix A) plus the leading two-loop thresholds 0{asQ:t) 

and Ola'f). The former is compnted for generic SUSY spectra while the latter (which is 
generically smaller) is only compnted for degenerate scalars. 

The hnal expression for the Higgs mass can thus be written as 

ml = v^[X{mt) + 6X{mt)], (1) 

where v = 246.22 GeV and SX{mt) are the SM threshold corrections (here compnted np to two 
loops) to match the Higgs pole mass to the MS rnnning qnartic conpling. The conpling X{mt) is 

^ As usual we interpreted the experimental value rut = 173.34±0.76 [31] as the pole mass, systematic uncertainties 
coming from this choice can be estimated by rescaling the experimental error on the top mass. 

^Since we do not perform a complete N^LO computation this last correction is also used to evaluate the 
uncertainties from higher order terms. 

^As for the three-loop top Yukawa threshold, the four-loop QCD corrections [33,34] to the strong coupling RGE 
has been used to estimate the uncertainties. 
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derived using the RGE and the boundary conditions at the SUSY scale (see below). The RGE for 
the Higgs quartic coupling are solved together with gauge and Yukawa couplings at three loops, 
in particular the top Yukawa yt is extracted from 


mt = + 5yt{mt)) , 


( 2 ) 


where 6yt{mt) is the SM threshold correction matching the top pole mass with the MS top Yukawa 
coupling and here computed at NNLO, and N^LO in the strong coupling. The matching at the 
SUSY scale Q is instead given by 

^ 9^(0+9'^W) + AA'‘> + AA>y. + AA® , (3) 

where contains the 1-loop thresholds matching the Higgs quartic coupling \{Q) in the MS- 

scheme with the one computed in full SUSY in terms of soft terms and couplings, in the DR or 
OS schemes. AA^^^^ are the leading two loop threshold corrections further discussed below. 

If some of the superpartners are light compared to the rest of the SUSY spectrum, as in the 
case of Split SUSY, a new mass threshold develops. In this case two matchings are in order, the 
hrst at the Split scale between the SM and the Split SUSY theory, and the second at the SUSY 
scale, between the Split theory and the MSSM. The evolution up to the Split scale is the same 
as in the previous case. We then used 1-loop thresholds to do the matchings and 2-loop RGEs 
to rnn the Split-SUSY theory [8,23]. We will show in section 3.4 that the simplest approach also 
works in the Split case, i.e. the effect coming from the splitting of the fermions from the scalar 
superpartners do not need RGE resnmmation in the parameter region relevant for the observed 
Higgs mass. 

Onr computation is very close to the one in [26], in particular we added the contributions 
from the bottom and tan Yukawas, relevant at large values of tan/3, we recomputed the two-loop 
thresholds 0{atas) using the effective potential in [35], and we also included 0{a^) corrections 
computed for degenerate scalar masses. 

The general expression for the two-loop 0{atas) corrections is too long to be reported here, 
but can be accessed through the computer code SusyHD provided in [27] for DR and OS schemes. 
Our computation in the DR scheme agrees^ with the one of [26]. In the limit = mu^ = mi 
and vanishing gluino mass the OS expression takes the simple form 


aa(2 ) = 
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while for M 3 = mq^ = m^/g = mi 


A 3(2) _ VtOl 
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( 5 ) 


'^We thank the authors of [26] for providing the explicit expression of their result for the cross-check. 
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where Xt = Xt/mi, Xt = At — 11 / tan 13 and Q is the renormalization scale. The dehnition we used 
for Xt in the OS scheme is given in eq. (27) in appendix A. 

The two-loop 0{at) supersymmetric threshold correction to the quartic coupling can be de¬ 
rived from the corresponding correction to the Higgs mass. We derived it under the simplifying 
assumption of degenerate scalars while the fi parameter and the renormalization scale are left free. 
We used the results in ref. [15] for the 0{a^) correction to the Higgs mass calculated using the 
effective potential technique in DR. Converting in the one-loop 0{at) correction the DR superpo¬ 
tential top Yukawa coupling and MSSM Higgs vev into the MS SM top Yukawa and EW vev will 
produce an additional shift contribution at two loops. Analogously for the OS computation there 
is an extra shift from converting the stop masses and mixings in the one-loop corrections. We 
subtract the 0{a^) top-quark contribution because it already appears in the matching at the EW 
scale. Finally, it is important to notice that there is also a contribution to the matching of the 
Higgs mass (and the quartic coupling) at the SUSY scale induced by the one-loop contribution of 
the stops to the wave-function renormalization of the Higgs held, which is instead absent in the 
0{atas) corrections. The complete expression with the details of the calculation can be found in 
appendix A, a simplihed expression in the OS scheme for the case ^ = mi and large tan [3 reads 


= 


(47r)^ 


3+fY 


—Xt + -X® 
6 * 6 * 


3 * 




As a cross check we verihed analytically that the two-loop 0{atas) and 0{a^) threshold corrections 
to the quartic coupling (under the assumption of degenerate scalars) cancel the dependence on 
the renormalization scale of the Higgs mass at the same order. 

Finally we also included the 1-loop threshold corrections from the bottom (and tan) sector, 
which are relevant in the large tan/3 region. The explicit expressions can be found in the ap¬ 
pendix A. At large tan /3, depending on the size and sign of other parameters, such as the fi term 
and the gaugino masses, the net effect is that of reducing the value of the Higgs mass. This effect 
may even cancel the tree-level contribution and allow for larger SUSY scales (see section 3.3). 

The relevance of the supersymmetric thresholds decreases as the SUSY scale increases because 
of the evolution of the SM running couplings. Among the missing SUSY threshold corrections 
the most important are the two-loop 0{at) when the scalars are not degenerate, the two-loop 
0{ata,asQ:) proportional to the electroweak gauge couplings and the three-loop 0{ata‘j.). In the 
case of large tan (3 and sizable fi parameter, the two-loop corrections proportional to the bottom 
Yukawa can also be relevant, they include the O^abds, abdt, oil) corrections. The contribution 
of the missing SUSY thresholds to the Higgs mass is estimated to be below 1 GeV even for a 
spectrum of superparticles as low as 1 TeV, see the next section. 


2.2 Estimate of the uncertainties 

In the EFT approach to the calculation of the Higgs mass in SUSY, the uncertainties can be recast 
into three different groups: 

1. SM uncertainties: from the missing higher order corrections in the matching of SM couplings 
at the EW scale and their RG evolution; 
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Figure 1: Breakdown of the uncertainties for a 125 GeV Higgs mass as a function of the (degenerate) superparticle 
masses tosusy- The Higgs mass has been kept fixed at 125 GeV by varying either the stop mixing (with fixed 
tan(5 = 20 for msusY < 20 TeV, left panel of the plot) or tan (5 (with vanishing stop mixing for msuSY > 20 TeV, 
right panel of the plot. Note that for tosusy <2 TeV (the gray region) the 125 GeV value for the Higgs mass cannot 
be reproduced anymore but is within the theoretical uncertainties. The black “total” line is the linear sum of the 
theoretical uncertainties from SM, SUSY and EFT corrections (in dashed lines). The dotted line A^p corresponds 
to the 2a experimental uncertainty on the top mass. 


2. SUSY uncertainties: from missing higher order corrections in the matching with the SUSY 
theory at the high scale; 

3. EFT uncertainties: from missing higher order corrections from higher dimensional operators 
in the SM EFT and other EW suppressed corrections C>(n^/mgusY)- 

Fig. 1 summarizes the importance of the individual sources of uncertainty as a function of the 
SUSY scale. For dehniteness we took the superpartners to be degenerate with mass msusY, the 
Higgs mass has been kept hxed at 125 GeV by varying either the stop mixing (with fixed tan /3 = 20 
for msusY < 20 TeV) or tan/3 (with vanishing stop mixing for msusY > 20 TeV). We will now 
discuss these uncertainties individually. 


SM uncertainties 

As described in the previous section, in our computation we employed full SM three-loop RGE 
and two-loop matching conditions at the EW scale to relate the pole masses my and mt and the 
gauge couplings to the MS running couplings at the high scale. We also included the 3-loop O(ttf) 
corrections to the top mass matching. This is expected to be the leading higher-order correction 
and the missing 3-loop matching and 4-loop running corrections are not expected to give larger 


7 









> 

0 ) 

S 

< 


-3 

0.5 


V 

\ m 

\2o 

N 

N 

N 

N 

0(ci:,as+a,h ^ v 

susY~^ TgV, 20 

X(= ^SUSY 

- 



Q / /Dsusy 



Figure 2: The uncanceled scale dependence from higher order corrections is largest at maximal stop mixing and 
small stop masses (here taken 1 TeV). Including only 1-loop SUSY threshold it amounts to up to a 3 GeV shift 
of the Higgs mass, when the scale is changed by a factor of 2. It reduces to below 1 GeV when the leading 2-loop 
0{atOts) and 0{a1) corrections are included. 


effects. Still, we conservatively used the 3-loop 0{a^) corrections to estimate the SM uncertainties 
from the higher-order missing corrections, although the latter are probably smaller^. 

The full SM uncertainty in fig. 1 has been computed by summing the effects from 0{al) 
corrections to the top mass and the ~0.15 GeV estimate [30] of the 3 loop corrections to m^. 
While the latter corrections would be formally of the same order as the corrections induced by 
the 2-loop SM corrections to the matching of the top Yukawa in a hxed order computation, in the 
EFT approach they are actually subleading because the 2-loop corrections from the top sector get 
RG enhanced. The net effect from these SM corrections amount to a shift to the Higgs mass of 
order 0.5-^l GeV for msusY ~ 1 10^ TeV. The uncertainty slowly increases with the SUSY scale 

as a result of the longer RGE running. 

SUSY uncertainties 

The matching between the SUSY soft parameters and the SM couplings includes full one-loop 
threshold corrections (including also bottom and tan Yukawa corrections) plus the leading two- 
loop corrections O(atas) and 0{a^) (the latter only in the simplihed case of degenerate scalar 
masses). While 0{atas) can give large effects, the corrections from Olaf) are substantially smaller 
and other missing 2-loop thresholds are expected to be even smaller. Since a missing threshold 
produce an uncanceled renormalization scale dependence in the final Higgs mass, such dependence 
can be used to estimate the missing corrections. In the worst case (maximal stop mixing and small 
SUSY scale) the uncanceled scale dependence from the 1-loop thresholds may shift the Higgs mass 
by roughly 3 GeV when the renormalization scale is changed by a factor of 2. 0{atas) reduce the 


^For this reason our theoretical uncertainty from the SM calculation is somewhat larger than the one quoted 
for example in [30], which uses the same precision for the computation of the stability of the SM Higgs potential. 





















shift® to 1 GeV and O^a^) further down below 1 GeV, see £g. 2. 

The uncertainty from SUSY thresholds in fig. 1 has been estimated by taking the maximum 
of the shifts induced by varying the SUSY matching scale by a factor of 2 or 1/2 with respect to 
'fTT'SVSY- The impact on the uncertainties from missing SUSY threshold corrections greatly reduces 
away from maximal stop mixing and when the stop masses are increased, the latter effects due 
to the reduction of the SM couplings from the RGE evolution. This fact is manifest in the EFT 
approach, less so in others, which require a careful resummation of logs. 


EFT uncertainties 

The last source of uncertainties is intrinsic to the EFT approach and comes from neglecting higher 
dimensional operators below the SUSY scale. Such corrections decouple fast, as powers of the EW 
scale over the SUSY scale, however they could become important for light SUSY scale. Given the 
relatively high value of the Higgs mass and the bounds from the LHG, superpartners are expected 
to lie above the TeV scale, reducing the relevance of these corrections only to very particular 
corners of the parameter space. 

At tree level the only source of power corrections comes from the heavy Higgs states—if rriA is 
close to the EW scale the mixing effects may become important and the tree-level expression in 
(3) gets corrections of order 

2 

6eft^ — —y sin^(2/3) -+-.... (7) 

m\ 

Such contributions become important only when mA is particularly light (m^ ^ 200 GeV), in 
a region of the parameter space which is already disfavored by indirect bounds on the Higgs 
couplings [36]. 

Gorrections from higher dimensional operators induced by the other superpartners enter only 
at one loop (such as the other scalars and the EWinos) or at two loops (gluino). The most 
dangerous corrections are thus expected to come from the stops, they are of and can 

get enhanced at large stop mixing. We estimated such corrections by multiplying the one-loop 
corrections by v'^ Numerically, for stops above 1 TeV, even at maximal mixing, these 
corrections are below 1 GeV and rapidly decouple for heavier stop masses, see fig 1. Having 
lighter stops may require to take such corrections into account, although they start being too 
light to accommodate the observed value of the Higgs mass (see sec. 3.2). Gonsistently with such 
approximation we also neglected EW corrections to the sparticle spectrum. 

The total EFT uncertainty in £g. 1 has been estimated by taking the sum of the single con¬ 
tributions to AA from each SUSY particle with mass mi multiplied by the corresponding factor 

In conclusion, for stops above the TeV scale power corrections are small, justifying the use of 
the EFT. 

®In the right plot of fig. 2 the scale dependence left after the inclusion of the 0{atas) corrections seems to be 
smaller because of an accidental cancellation for those particular values of the parameters. 
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Combined uncertainties 


Fig. 1 summarizes the impact of the various uncertainties to the determination of the Higgs mass 
as a function of the SUSY scale, tan/? and the stop mixing, in the relevant region of parameters 
that reproduces the measured value of the Higgs mass. For definiteness we took a degenerate 
spectrum of superpartners, we checked that the size of the uncertainties remains of the same 
order when this assumption is relaxed. The dominant source of error comes from higher order 
corrections in the matching and running of the SM couplings. SUSY thresholds are only important 
for low SUSY scale and large stop mixing, while power corrections are negligible throughout the 
parameter space unless some of the sparticles are very close to the EW scale. 

It is fair to say that, for most part of the relevant parameter space, the Higgs mass in the 
MSSM has reached the same level of accuracy as the determination of the Higgs potential in the 
SM. Further improvements from the theory side can be achieved by extending the SM calculations 
at higher orders. The size of the uncertainties remains practically unchanged in the split scenario, 
where the fermions are parametrically lighter than the scalar superpartners. 

The total theoretical uncertainty (computed here conservatively as the linear sum of the three 
sources of errors discussed above) is of order 1 GeV or below for most of the parameter space. 
It is thus below the error induced by the experimental uncertainty in the value of the top mass. 
Indeed, the latter produces a shift in the Higgs mass of order 1.5-^2.5 GeV depending on msusv, 
when the top mass value is changed by 2a=1.5 GeV. The error increases with msusv due to RGE 
effects. 

As usual, estimates of theoretical errors provide only for the order of magnitude of the expected 
corrections and must be taken with a grain of salt. However since for most of the parameter space 
the error is dominated by the SM uncertainties, where we have been rather conservative, the 
estimate of £g. 1 should represent a fair assessment, at least away from the lower end. 

2.3 Comparison with existing computations 

Our EFT computation agrees within the uncertainties with all the others which use the same 
technique. As already noticed in [26], however, the EFT computation seems to give a smaller 
Higgs mass with respect to other approaches, such as those based on full diagrammatic and 
effective potential computations such as [37-40]. In some cases the disagreement amounts to up 
to ~10 GeV, well beyond the expected quoted uncertainties, even in regions of parameter space 
where both approaches are expected to hold. 

A comparison between the EFT computation and some of the available computer codes is 
shown in fig. 3. The disagreement is around 3 GeV for msusv > TeV at large tan/? and zero stop 
mixing and increases up to 9 GeV for maximal mixing and msusv = 2 TeV. 

The large disagreement with the FeynHiggs 2.10.1 code can mostly be understood as follows. 
The computation in [41] included full 1-loop plus the leading 2-loop SUSY corrections of the 
Higgs mass with partial 2-loop RGE improvements. Gonsistently with this however, they did not 
include 2-loop corrections to the matching of the top Yukawa coupling. Instead, the use of the 
N^LO formula shifts the top mass by roughly 4 GeV. Hence, the bulk of the disagreement seems 
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Figure 3: Comparison between the EFT computation (lower blue band) and two existing codes: FeynHiggs [ 4 I] 
and Suspect [39]. We used a degenerate SUSY spectrum with mass tosusy in the DH-scheme with tan/3 = 20. 
The plot on the left is mh vs msvSY for vanishing stop mixing. The plot on the right is mn vs Xt/msuSY for 
uisuSY = 2 TeV. On the left plot the instability of the non-EFT codes at large tosusy is visible. 


due to the missing 2-loop corrections in the top mass^. Note that, as discussed in the previous 
section, the uncertainty in the EFT approach is dominated by the 3-loop top matching conditions, 
the 2-loop ones are thus mandatory in any precision computation of the Higgs mass. We checked 
that after their inclusion, the FeynHiggs code would perfectly agree with the EFT computation 
at zero squark mixing. At maximal mixing the disagreement would be reduced to 4 GeV, which 
should be within the expected theoretical uncertainties of the diagrammatic computation. 

For comparison, in Eg. 3 we also show the results obtained with a different code (Suspect [39]) 
which uses a diagrammatic approach but unlike FeynHiggs, does not perform RGE improvement 
and its applicability becomes questionable for msusY in the multi TeV region. 


3 Results 

After having seen that the EFT computation is reliable for most of the relevant parameter space 
we present here some of the implications for the supersymmetric spectrum. Given the generic 
agreement with previous computations using the same approach, we tried to be as complemen¬ 
tary as possible in the presentation of our results, putting emphasis on the improvements of our 
computation and novel analysis in the EFT approach. 

3.1 Where is SUSY? 

Fig. 4 represents the parameter space compatible with the experimental value of the Higgs mass in 
the plane of {mi/ 2 , mo) for zero (blue) and increasing values (red) of the stop mixing. For simplicity 
we took degenerate scalar masses mo as well as degenerate fermion masses mi /2 = Vfi 2,3 = ft. All 

^It was brought to our attention that a similar observation was also made in [42]. 
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Figure 4: Higgs mass constraint on the value of scalar (mo) and fermionic (mi/ 2 ) superpartners (taken de¬ 
generate). The upper blue region refer to zero mixing, the lower red to increasing values of the stop mixing. The 
lighter bands corresponds to the uncertainty from the top mass. The gray shaded region corresponds to tachyonic 
on-shell masses for the squarks. The non-vertical thin dashed lines correspond to on-shell values for the squark 
masses: when toi /2 grows the DR mass mo must be increased to keep the on-shell mass constant. The thick dashed 
burgundy lines correspond roughly to the expected reach of LHCIA and of an hypothetical 100 TeV machine. 


SUSY parameters of this plot are in the DR scheme®. The figure highlights a number of features: 

• The main effect at small fermion masses is given by the scale of the scalars (in particular 
the stops). The lower part of the allowed region corresponds to large values of tan/3 > 10. 
Lowering tan/3 allows to access larger scalar masses (see also fig. 6 below). 

• The dependence on the fermion masses can be understood as follows. For mi /2 ^ trtQ the 
biggest contribution comes from the higgsino-wino loop in the running of the Higgs quartic. 
It makes the quartic coupling run larger in the IR thus making the Higgs heavier. This 
correction is only there when both wino and higgsino become light. There is also a smaller 
1-loop correction from the individual EWinos, which affects the running of the EW gauge 
couplings. They make the gauge coupling run bigger in the UV increasing the tree-level 

®A11 DR parameters are computed at the scale Q = unless specified otherwise. 
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contribution to the Higgs quartic (3) and thus its pole mass. Lowering the gluino mass 
decreases the Higgs mass but the effect is two-loop suppressed and only non-negligible at 
large stop mixings. The region mi /2 ^ 2mo should be treated with care. In the DR scheme 
there are negative quadratic corrections to the squark masses proportional to the gaugino 
masses [43] 




, Aas , ^ 


log I ^ I - 1 


+ ... 


( 8 ) 


In particular when M 3 becomes larger than roughly a factor of four with respect to the squark 
masses the corresponding on-shell masses become tachyonic. Just before this happens the on- 
shell masses (the dashed lines in the hgure) start becoming smaller and smaller with respect 
to the DR parameters, in this tuned region large corrections make the DR computation 
unstable. This explains the strong apparent dependence on mi /2 on the right-hand part of 
the plot, which would disappear if plotted in terms of the on-shell masses. We decided to 
keep the plot in terms of the DR parameters to highlight the tuning required to explore such 
region. 


• Current LHC searches already probed squark and gluino masses up to 1.5 TeV circa [44]. 
This corresponds to the very lowest part of the allowed parameter space, where the stop 
mixing is maximal, tan is large and fermions must be lighter than scalars. This, of course, 
with the caveat that the strongest experimental bounds apply to hrst generation squarks 
and gluino while the Higgs mass mostly depend on the stops and (somewhat weaklier) on 
EWino. With the same caveat LHC14 should eventually be able to more conhdently explore 
the same region (extending the squark-gluino reach to 3 TeV, see e.g. [45]), while the small 
stop mixing region could only be reached directly with a 100 TeV machine (capable of 
probing colored sparticles of roughly 15 TeV masses, see e.g. [45]). Of course (mini-)Split 
scenarios where the heavy scalars are responsible for the Higgs mass and the light fermions 
are within reach at lower energies remain a valid possibility. 


3.2 The EFT gets on-shell 

Previous computations using the EFT approach have used the DR scheme for the SUSY and 
the soft parameters. This scheme has the advantage of being the natural framework for the 
computations of the soft parameters in theories of SUSY breaking. In some cases, however, it 
results inadequate for the computation of the Higgs mass. 

First of all, physical on-shell masses are needed to compare theoretical computation with 
experiments. While the difference in the schemes is one-loop suppressed, there are non-decoupling 
effects which require care. For example the difference between the on-shell and the DR squared 
mass of the squarks receives an additive one-loop correction proportional to the gluino mass 
squared, see eq. ( 8 ). Such correction is negative and big—it is enough for the gluino mass to be 
a factor of four above the squark masses to drive the corresponding on-shell mass tachyonic. 

For similar reasons in the DR scheme the gluino contribution to the Higgs mass does not 
decouple [16]. Another consequence is the instability of the Higgs mass with respect to the 
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Figure 5: Allowed values of the OS stop mass reproducing ruh = 125 GeV as a function of the stop mixing, with 
tan j3 = 20, p, = 300 GeV and all the other sparticles at 2 TeV. The hand reproduce the theoretical uncertainties 
while the dashed line the 2a experimental uncertainty from the top mass. The wiggle around the positive maximal 
mixing point is due to the physical threshold when crosses + mt ■ 


renormalization scale—even if the on-shell squark masses are positive, the DR stop mass becomes 
highly sensitive to the renormalization scale when the gluino is more than a factor of 2-^3 above 
it, which results in an instability of the estimate of the Higgs mass. What is happening is that 
the physical on-shell squark masses becomes tuned and highly sensitive to the soft parameters. 
The situation is similar to trying to compute the Higgs mass in terms of the soft parameter rn\j^ 
instead of the on-shell (tuned) EW vev v. 

All these problems disappear in the OS scheme, the gluino decouples up to a physical log 
correction [16], there are no tachyons since the physical OS masses are given as input and larger 
hierarchies can be introduced safely within the SUSY spectrum (with the usual caveat that large 
logarithms may require resummation). Besides, the input masses are directly the physical quan¬ 
tities to be compared with experiments. 

For these reasons we also performed our computation in the OS scheme. Fig. 5 shows an 
application of such calculation. It corresponds to the region of allowed OS stop masses (taken 
degenerate in this case) which reproduces the observed Higgs mass for different A^-terms. Our 
dehnition of At in the on-shell scheme, eq. (27), is different from the usual one, this explains why 
the point of maximal mixing is not at Xt/mi ~ 2. In the spirit of natural SUSY [46-48] we kept 
the higgsino light at 300 GeV while the gauginos and hrst generation squarks safely above collider 
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bounds at 2 TeV. The lightest stop masses allowed in this case (for maximal stop mixing) are 
about 1.7±0.4 TeV, in the region where the EFT approach should be reliable. 

Had we drawn the same plot in terms of the DR masses we would not be able to draw the 
same conclusion—the error would blow up in the region where the stops are sufficiently lighter 
than the gluino. 


3.3 Large-tan ,5 High-Scale SUSY strikes back? 


While for most values of tan [3 the contributions from the bottom and tau sector can be neglected, 
at very large tan/? the corresponding superpotential couplings become large and their effects to 
the SUSY threshold can eventually dominate over the others. In particular the one loop sbottom 
threshold to the Higgs quartic coupling at leading order in tan f3 and degenerate sbottoms (mg) 
reads 

V,4 ,,4 

(9) 


AAy = 

b 


yt 


y 


327r2 mf ’ 

b 


and analogously for the tau. At tree-level the superpotential Yukawa coupling yi, is related to the 
SM Yukawa yb by jjb = yb/ cos 13. At large tan f3, jjb may become larger than one. In this situation 
the negative threshold correction (9) may cancel or even overcome the tree-level contribution, 
especially at high SUSY scale where the SM EW gauge couplings are smaller. This effect may 
allow to maintain rrih = 125 GeV and large tan/3 with arbitrary heavy scalar fields, reopening 
the High Scale SUSY window above 10^° GeV, which was thought to be excluded by the Higgs 
mass within the MSSM. As an example, we show in Eg. 6 how the msusY-'vs-tan /3 plot would 
look like at large tan (3 after including the leading bottom (and tau) contributions. For the plot 
we chose degenerate spectrum with mass msusv, negative y = —msusv and At = ?71susy/2. 
The SUSY parameters are given in the OS scheme. The behavior at small and moderately large 
tan f3 (tan f3 < 40) is well-known [23,26,28]. However further increasing tan [3, the bottom coupling 
dfe = yl/{‘in) grows, decreasing the Higgs mass [49]. For very large tan [3 the tree-level contribution 
to the bottom mass is so snppressed that loop corrections cannot be neglected [50,51]. In fact the 
bottom mass receives corrections from SUSY breaking proportional to Vu = nsin/?, i.e. not tan [3 
snppressed 


yb = yb cos /3 + yb sin /3 


8 a, /4M3 W 1 

3 47r m? v m? / in sin^ [3 rri^ \mt / 

b b ' b b 


+ ... 


( 10 ) 


where F{x) = (1 — a; -|- a;loga;)/(l — x)^ and we considered uiq^ = mus = niDs = The 

loop corrections are proportional to y and a combination of ganginos and A-terms. If the latter 
are small or have opposite sign with respect to y (like in Eg. 6), ijb will become strong at large 
tan(3 in order to reprodnce the observed bottom Yukawa (the red region in the plot). Before 
that, ijb is large enongh to make the threshold (9) win over the tree-level contribntion and allow 
^susY > 10^° GeV at large tan (3. For example the observed Higgs mass can be reproduced for 
GUT scale SUSY with tan(3 ~ 200. In this case the bottom conpling ab ~ 0.5, which is still 
pertnrbative bnt with a very close Landan pole A^p ^ lOmsusY- The pertnrbativity of Ab could 
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Figure 6: Effects of the sbottom threshold at large tan/3. The blue band corresponds to the ruh = 125 GeV 
constraint (the width is given by the estimated theoretical uncertainties) for different values of tan /3 and the degen¬ 
erate SUSY mass rnguSY - IFe fixed p, = — tosusy o-nd At = msusY/2. Dashed and dot dashed lines correspond to 
different values of the bottom and tau couplings respectively. The red region correspond to tachyonic Higgs and/or 
non-perturbative bottom Yukawa coupling. 


be improved by choosing larger /i terms, however this may become in tension with bounds from 
tunneling into charge/color breaking vacua [52,53]. We do not know what are the corresponding 
bounds on the fx term in this regime, this require a dedicated study which is beyond the scope of 
this work. 

We thus hnd that the upper bound of 10^° GeV on the SUSY scale from the observed Higgs 
mass may not apply for arbitrary values of tan/3 but only for small to moderately large tan/3. 
High scale SUSY at larger tan/3, however, requires large /i terms, gauginos may be lighter but 
not too much since they receive loop corrections. Therefore high scale Split SUSY does not seem 
possible in this way. 

If gaugino masses and/or H-terms are large and with the same sign as p, the loop corrections 
may saturate the full contribution to the physical fermion mass. If this happens, arbitrary large 
values of tan /3 can be reached without ever running into strong coupling effects. 

Finally for smaller p (not shown in the plot) the bottom-tan sector remains decoupled from the 
low energy Higgs, the threshold (9) is never important and msusY at large tan/3 stays constant. 
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Figure 7: Constraints in the (mQ,mi/ 2 ) plane from nih = 125 ± 1 GeV^ for different values on tan/3 using 

the proper split-SUSY computation (dark yellow, dashed line) or the approximate high-scale computation (red, 
continuous line) which does not resum the logs induced by the splitting of mi /2 from mg. The agreement is 
remarkable in the whole relevant parameter space. 


The bottom and tau Yukawa couplings still become strongly coupled at large tan (3 but the effect 
on the Higgs mass remain small. Of course the effect from the new physics present at the strong 
coupling scale is model dependent and may be important. 

3.4 Split vs High-Scale SUSY computation 

As mentioned before, in (mini)split-SUSY scenarios, where gauginos and possibly the higgsinos 
are sensibly lighter than the scalar sector, a new mass scale is present and large logarithms 
may require resummation. In this case the correct procedure would be: 1) to interrupt the SM 
running at the split scale, where the light fermion superpartners are, 2) to match to the split- 
SUSY effective theory, which includes SM particles and the fermion superpartners, 3) to perform 
a second running within the new EFT and eventually 4) to match to the full SUSY theory at 
the scalar mass scale. This procedure, which has been employed since the birth of split SUSY, 
and became more popular recently after the Higgs discovery, is more involved than the high-scale 
SUSY computation. Besides, the thresholds and the RGE of the split EFT are only known at a 
lower order in perturbation theory. 
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Note however the following. The leading effect of resumming the logarithms generated from 
splitting the fermion superpartners is to change the running of the Higgs quartic and EW gauge 
couplings. Numerically the change in the running of the Higgs quartic coupling is the leading 
contribution but it is only present when both higgsino and gauginos are light. The effect from the 
change in the RGE of the EW gauge couplings is instead smaller and it is further suppressed at 
small tan/3, exactly when the logarithms are the largest. 

The observed value of the Higgs mass is not very large and its value limits the SUSY scale to 
roughly 10^° GeV. This scale gets further reduced to GeV if the SUSY fermions are split, 

as an effect of the extra contribution to the running of the Higgs quartic. This translates into an 
upper bound on how large the logarithmic thresholds from splitting the fermions can grow. 

It turns out that in the whole parameter space relevant for the observed value of the Higgs 
mass, the effect of resumming the logs of the splitting between fermion and scalar superpartners is 
negligible and the results obtained with the single-scale SUSY theory are reliable. This is shown 
in fig. 7 where we compare the computation made resumming the logs of the split threshold, using 
an intermediate Split-SUSY EFT, with the one that does not resum the logs, which uses one scale 
only and just the SM RGE up to the scalar masses. The agreement between the two procedures 
is impressive. In the worst case (tan/3 = 1, fermions at 200 GeV and scalars at 10® GeV) the 
mismatch is less than 1 GeV, well within the estimated uncertainties. It can also be seen that 
the two procedures start deviating exactly at that point. Indeed, had the SUSY scale and the 
splitting between fermions and scalars been bigger, the two computations would start deviating 
sensibly, fortunately that region is not relevant for ruh = 125 GeV. 

We conclude that for all the relevant parameter space the computation of the Higgs mass can 
reliably be made using only the SM as EFT up to the scalar masses, independently of the scale 
of the fermions®, whose main effect is well approximated by the one-loop thresholds at the SUSY 
scale. 


4 The SusyHD code 

The computation described in the previous sections has been implemented into a simple Math- 
ematica [56] package, SusyHD (SUperSYmmetric Higgs mass Determination), which we made 
public [27]. The package provides two main functions that compute the Higgs mass and its theo¬ 
retical uncertainties from the input soft parameters, and an auxiliary function to change the SM 
parameters {rrit and q;^). 

The most time consuming part of the EFT calculation is the integration of the RGE. The 
code avoids such step by using an interpolating formula for the solution of the RGE, which is 
only function of the amount of running \og{Q/mt) and the value of the Higgs quartic coupling at 
the high scale set by the SUSY threshold corrections, eq. (3). The interpolating formula 

only depends on the SM parameters, so the RGE integration needs to be run only once, when 
the package is hrst called or if the SM parameters are changed. The result is a very fast code 

®The usual caveat from u/msusY corrections applies when the fermions are very close to the EW scale. In this 
case the full contributions from the SM-I-fermion states [54,55] should be used in the matching at the low scale. 
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which allows to effectively use the observed value of the Higgs mass as a constraint for the SUSY 
parameter space. All plots of this paper have been generated with SusyHD. 

The input SUSY parameters can be given in either DR or OS schemes and thanks to the EFT 
approach they can be arbitrary heavy. The code also accepts simplified input where not all the 
SUSY parameters needs to be specified. There are also extra options which allow: 1) to switch 
off independently some of the higher order corrections, 2) to change the matching scale Q, 3) to 
use the full numerical code, which integrates the RGE numerically and 4) to use the Split SUSY 
code which integrates the RGE in two steps: SM up to the fermion scale and Split-SUSY up to 
the scalars. The function that computes the theoretical uncertainties accepts also the option to 
compute the individual uncertainties coming from the SM corrections, the SUSY thresholds, and 
the EFT approximation. 

All the necessary documentation can be downloaded with the code from [27]. 


5 Phenomenological Applications 

5.1 Predicting the spectrum of Minimal Gauge Mediation 

Gauge mediated supersymmetry breaking (GMSB) [57-64] is among the simplest and most elegant 
calculable mechanisms for generating the MSSM soft terms. A very special property is the absence 
of dangerous FGNG, a very rare property in extensions of the SM Higgs sector, supersymmetric 
and non. 

However, using GMSB to implement a natural solution of the hierarchy problem has always 
been hard. The main obstruction being the /i problem, viz. why the supersymmetric higgsino 
mass happens to be at the same scale of the SUSY-breaking soft terms. Solutions of the fi problem 
generically produce a /i/R/i problem [65]: both /i and Bfi are generated radiatively at the same 
order in perturbation theory, which produces an unwanted hierarchy, Bfi ^ /i^. Solutions to the 
/i/R/i problem exist (see e.g. [64-66] and references therein) but at the cost of an excessive model 
building. 

All these problems arise when we try to obtain a natural SUSY spectrum—it is like we are 
forcing the theory to do something it was not meant to. In line with what discussed in the 
introduction, we are then going to relax this requirement and try to use experiments instead of 
naturalness to infer the properties of new physics. 

The apparent gap between the EW and the new physics scale motivates us to revisit the 
simplest and more elegant GMSB model, minimal gauge mediation (MGM)^°, without the unnec¬ 
essary baroque model building associated to the Higgs sector. Indeed, ignoring the naturalness 
problem allows us to also ignore the /i-problem, as the two are closely related. When /i is much 
larger than the soft masses EWSB is not possible, when fi is much smaller, the EW scale v would 
be of order the soft masses. Therefore if the SUSY scale is above v, fi must automatically be close 
to the SUSY scale in order for EWSB to be tuned to its experimental value. 

^°Here by MGM we really mean the most minimal realization, where the Higgs sector only receives the standard 
gauge mediated contribution, fi is & free parameter and is generated radiatively in the IR. 
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Figure 8: Prediction for the spectrum of MGM after imposing the constraint from the Higgs mass (or better from 
the top mass). For each superpartner we plot the allowed range of masses (in TeV) for four different combinations of 
-/V = 1(3) and M = 10^(10^^) TeV. For each mass the lowest (highest) value corresponds to increasing (decreasing) 
the value of the top mass by 2a with respect to its experimental central value. The values of tan /3 at the bottom 
(top) side of each of the four bands, from left to right, are 58 (42), 49 (45), 56 (29) and 44 (45) respectively. The 
three differently shaded areas represent “pictorially” the existing LHC8 bounds and the expected reach at LHClf 
and at a future 100 TeV collider, respectively from the bottom. 


In MGM all soft masses are generated with the same order of magnitude by the gauge mediated 
contribution, one gauge loop below the scale A = F/M (the ratio between the effective scale of 
SUSY breaking F and the mass of the messengers). Besides A, the spectrum also depends, in a 
milder way, on the actual mass of the messengers M, which determines the amount of running 
of the soft parameters, and the number of messengers N (typically Y = 1 or 3 for a vector like 
messenger in the 5 or 10 of SU(5) respectively). 

As mentioned before, the /i-term, being supersymmetric, would be an independent parameter, 
but its value is hxed by requiring (tuning) the correct EWSB. Finally the A-terms and are 
generated radiatively from RGE effects. This fact has very interesting consequences [67,68]. First, 
being A and Bjx terms generated at the quantum level from gaugino masses and /x-term implies 
that the corresponding GP phases vanish, avoiding potentially dangerous bounds from EDMs. 
Second, small suppressed A-terms imply that the stop mixing will never be large, while small Bfx 
implies large values of tan (f. These two predictions combined with the measured value of the Higgs 
mass allows to £x also the overall scale A, which must then lie at around the PeV scale to produce 
the (P(IO) TeV SUSY scale required by the Higgs mass. The only remaining free parameters are 
the messenger mass scale M and their number N, which affect the properties of the spectrum in 
a milder way. 
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Using our computation for the Higgs mass we can thus predict the spectrum of MGM in 
terms of N and M, the result is shown in hg. 8. Four different spectra are reported, changing 
independently N (1 or 3) and the messenger scale M from M = 10^ GeV (to allow the use of 
leading 0{F/M) formulae) to M = 10^^ GeV (to avoid dangerous FGNG contributions from 
gravity mediated contributions). For each choice of N and M the spectrum is not completely 
determined because of the uncertainty in the Higgs mass computation. Indeed the effect of varying 
N and M is actually subleading with respect to the Higgs mass uncertainty. In the relevant region 
of parameters (msusY ~ 10^ TeV and small stop mixing) the Higgs mass determination is at its 
best (see £g. 1). Theoretical uncertainties are completely dominated by the SM ones, which are 
sub leading with respect to the experimental uncertainties in the top mass. In fact, what limits the 
prediction of the MGM spectrum is not the Higgs mass, or its determination in SUSY, but our 
poor knowledge of the top mass! Improvements in this quantity are required to further improve 
the predictions of hg. 8. The lowest (upper) bounds correspond to values of the top mass 2a 
above (below) its measured central value. The overall scale A results to be at the PeV scale, in 
particular it varies roughly from 0.5 to 2.6 PeV for different choices of the top mass, N and M. 
The values of tan f3 are typically around 45 but they can vary up to 60 and down to 30 in the 
corners of the parameter space, the corresponding values for the supersymmetric bottom and tan 
Yukawa couplings are largish (typically around 0.5-0.7) but remain always below the one of the 
top Yukawa. Similarly the stop mixing parameter is always small < 1. 

Except for the overall scale A, which is one order of magnitude larger than the one usually 
considered in the literature, the rest of the spectrum has the typical GMSB form, with bino or 
right-handed stau being the NLSP depending if Y = 1 or 3 respectively. 

On the experimental side, besides the simplified model and the generic SUSY searches, ATLAS 
and GMS also performed a number of dedicated GMSB searches [69-76], which exploits some of 
the most peculiar properties of its spectrum, such as photon and tans in the final states. Of 
particular relevance for this scenario is the direct search for the pseudoscalar Higgs boson A°, 
which, for the large values of tan/? predicted here, bounds mA ^ 800 GeV [76]. This channel 
appears to be the most powerful for MGM, with a slightly better reach than the standard GMSB 
candles. 

While a dedicated study is required, in fig. 8 we also show “pictorially” the existing experi¬ 
mental bounds and the expected reach at LHG14 and at an hypothetical 100 TeV machine. The 
latter are obtained by rescaling the pdf on the existing bounds [77, 78] and should serve only to 
guide the eye. However, given the expected scale of the spectrum we can confidently say that 
this model is mostly out of the reach of existing collider machines^^, but could be seriously (if not 
completely) explored by a 100 TeV hadron collider. In fact, MGM may well represent one of the 
strongest motivations for such machine. 

We checked that, while the values of tan [3 in this model are large, bounds from the rare decays 
Bd,s —t /i/i [79] are not strong enough to be sensitive to the spectrum in fig. 8 yet. An improvement 
on the experimental bounds by a factor 3-i-5 could be enough to start probing the bottom part of 
the spectrum. 

would like to point out that in some corners of the allowed parameter space A° may be light enough to be 
within the reach of LHC 14. 
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Figure 9: Allowed parameter space of minimal anomaly mediation in the plane (mo, p.) for different values of 
tan/?, after imposing the constraint from the Higgs mass. According to AMSB gauginos are one-loop lighter than 
the scalars, here taken with a common mass mo- The wiggle for negative p, at small tan/? is due to a cancellation 
in the one loop threshold correction from EWinos when fi crosses the gaugino masses. Values of tan fd > 3 are 
excluded by LHC bounds on Winos. The horizontal red line corresponds to M 2 — 3 TeV. 


In conclusion, MGM represents probably the simplest and most predictive implementation 
of SUSY. The whole spectrum is almost completely determined just by experimental data. In 
particular, the upper bound on the scale of the superpartners exists independently from any 
naturalness consideration, in fact the value of the Higgs mass predicts a SUSY scale not too far 
beyond our current reach. The model also makes other successful predictions such as: gauge 
coupling unihcation, the absence of SUSY particles at current hadron colliders, no EDMs and no 
deviation from any flavor observables. Gravitino may be dark matter although this possibility is 
more model dependent. 

5.2 Anomaly Mediation 

Minimal anomaly mediated supersymmetry breaking (AMSB) models [80, 81] probably provide 
for the simplest implementation of Mini Split supersymmetry. Scalars get their mass from gravity 
mediation of order the gravitino mass, while fermions, protected by R-symmetry, get one-loop 
suppressed soft masses from anomalies. The generation of the fi and Bfi parameters require also 
the breaking of the PQ symmetry so that the higgsino mass is practically a free parameter—it can 
be of order the gravitino mass or naturally smaller if the PQ breaking is not efficient. The theory 
is thus defined by three main parameters: the gravitino mass m 3/2 setting the scale of scalars and 
gauginos, the higgsino mass /x and tan fd which is determined by the details of the scalar masses 
and B^. 

Unlike in MGM the details of the scalar spectrum are model dependent, however, given the 
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large scales involved in this scenario, threshold corrections at the SUSY scale are almost irrelevant, 
for dehniteness we £x all the scalars degenerate, mo = m 3 / 2 . The actual value of the Higgs mass 
gives a further constraint on these parameters. It can be used, for example, to fix tan /3 in terms 
of the other two parameters. It is trivial to impose such constraint using SusyHD, the result is 
shown in fig 9. Values of tan /3 larger than 3-i-4 are already excluded, for they would require too low 
SUSY scale and the wino would lie below the LHC bounds [82]. Also a wino with mass ~ 3 TeV, 
which would provide for a good thermal dark matter candidate, corresponds to tan/? ~ 2 3. 

Bounds on such parameters from direct detection experiment can be found for instance in [78]. 

Since a LSP wino above ~ 3 TeV would overdose the universe in the minimal AMSB scenario, 
the allowed parameter space reduces to the “narrow” strip below the red line in fig. 9, and the 
one with \n\ < 1 TeV < M 3 where the higgsino is the LSP. Most of this parameter space could 
in principle be probed at a larger hadron collider and future dark matter experiments [78,83-87]. 
In this interesting region, tan/? is constrained between 2 and 3; scalars are clearly out of reach, 
between 10^ and 10^ TeV, but not heavy enough to guarantee the absence of FCNC [ 88 ]. Gauge 
coupling unihcation further prefers values of /i below O{10) TeV [28]. 


6 Conclusions 

We presented a calculation of the Higgs mass in the MSSM using the EFT approach, which 
improves previous computations by including extra two-loop SUSY threshold corrections, the 
contributions from the sbottom/stau sectors relevant at large tan/? and the implementation of the 
OS scheme, the relevant formulae can be found in the appendix A and in [27]. 

We also performed a study of the theoretical uncertainties, showing that for most of the relevant 
parameter space the error is sub-GeV and dominated by higher order SM corrections. The result 
is summarized in fig. 1 . 

The computation has been arranged into an efhcient computer package which we made publicly 
available [27]. The code exploits the power of the EFT approach, allowing to compute the Higgs 
mass for arbitrary heavy sparticles, even when a large hierarchy between fermions and scalars is 
present. Analytic formulae for the solution of the RGE make the code very fast, which allows to 
efficiently use the Higgs mass as a constraint on the spectrum. 

We then performed several studies on the implication of the Higgs mass constraint on SUSY; 

• In agreement with previous EFT computation we find that the SUSY spectrum needs to be 
a little heavier than expected, in particular stops below 2 TeV are disfavored (see fig. 5). 

• The upper bound on the SUSY spectrum, which is 0(10^°) GeV (0(10"^) GeV at large 
tau/d), can actually be relaxed without adding new degrees of freedom. At very large tan/?, 
if is not suppressed with respect to the scalar masses, sbottom/stau contributions may 
reduce the Higgs mass, allowing larger values for the SUSY scale (see fig. 6 ). 

• In mini-split SUSY, in the region of parameter space relevant for the Higgs mass, the effect 
of the thresholds from splitting the fermions from the scalars is completely captured by 
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the leading fixed order one loop corrections (see fig. 7). This allows to use the SM as an 
effective ffeld theory all the way up to the scalar mass scale, avoiding the need of using an 
intermediate split SUSY effective theory. 

• We point out that the value of the Higgs mass may be used to predict the spectrum of minimal 
gauge mediation, the simplest calculable SUSY model, almost completely. The spectrum of 
SUSY in this case can thus be bounded just by experimental data alone without the need 
of arguments based on naturalness. Interestingly enough the spectrum lies just above the 
expected reach of LHC14 (see fig. 8 ), making it an ideal target for a future 100 TeV hadron 
machine. 

• Finally we discuss about the analogous implications for anomaly mediation models, con¬ 
straining the allowed values of tan(3 and the scale of SUSY (see fig. 9). 

For most of the allowed parameter space the Higgs mass computation is dominated by the 
experimental uncertainty in the top mass. The theoretical uncertainties instead are mostly dom¬ 
inated by the SM higher order corrections. Only for maximal stop mixing and at the lightest 
possible stop masses uncertainties from SUSY corrections and from higher-order terms in the 
EFT expansion may become important. Improvements in this region can be achieved by includ¬ 
ing subleading two-loop threshold corrections neglected in this work, such as 0{aas^t) or O^ata"^), 
and the leading /rrigu^Y) corrections in the EFT expansion. 
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A SUSY thresholds 

This appendix is dedicated to some analytical expressions of the threshold corrections from inte¬ 
grating out supersymmetric particles that were not written in the body for the sake of read¬ 
ability. We start summarizing our conventions. For the numerical part we used the values 
rrit = 173.34 GeV [31], as{mz) = 0.1185 [89], ybim-t) = 0.0156 and yrirrit) = 0.0100 [25]. The 
MSSM Lagrangian is written with all the parameters in the DR scheme (or the “OS” scheme 
described below), and is matched with the SM Lagrangian with all couplings and masses in the 
MS scheme. For tan/? we used the definition of [26]. As in the rest of the paper, unless speci¬ 
fied otherwise, all the formulae are written in terms of the SM couplings (( 71 , 2 , 3 , yt,b,r and A, or 
ttj = gf/^Air) and at,b,r = ytbr/i^'^)) scheme and the soft parameters (masses and 

trilinear couplings) in the DR or OS schemes. 

The SUSY-breaking masses for the scalars of the i-th generation are denoted by uiq,, mjj., 
^Di, ^Li and niEi, the soft SUSY-breaking Higgs-squarks cubic couplings are written in terms of 
the superpotential Yukawas yt,b,T as at = ctb = {jb^b, ab = yb^b for the stops, sbottoms and 
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staus respectively, while the relative signs of the fi parameter, gaugino masses and A-terms are the 
same as in [90], so that the scalar mass mixings depend on Xt = At — /icot (3, = At, — /rtan/3 

and Xt- = A^- — fi tan (3. 

In this work, we extended the one-loop threshold in eq. (10) of [26] to inclnde also the 
tan /3—enhanced contribntions from integrating ont sbottoms and staus: 
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In the last three lines of the equation above we sum over the contributions of the stops, sbottoms 
and staus, where is the third component of weak isospin of the left-handed chiral multiplet 
to which the sfermions belongs, is the electric charge, X^ = {X^p/{mq^mu^), X^/{mq^mcs), 
X^KruL^mE^)}, Nf = {3,3,1} is the color factor, x^ = {rngg/m^/g,rngg/m^g,m^g/mEg}, and 
r-tj, = (1, ybCos/3/yb, Or cos(3/yr}. The latter coefficients take into account the tan/3 enhanced 
corrections discussed in sec. 3.3 which require resummation, the explicit expressions can be found 
e.g. in [25,27]. The loop functions are dehned in appendix A of [26]. Because of the smallness 
of the bottom Yukawa coupling, the one-loop 0{ab) SUSY threshold corrections are only sizable 
for large tan/3 and \y\ > yThg^rnb}. 

We obtained the two-loop 0{a() SUSY threshold corrections to the quartic coupling of the 
Higgs from the corresponding correction to the Higgs mass, under the simplifying assumption 
of degenerate scalars {mq^ = mu^ = rriA = 'mi) while the y parameter and the renormalization 
scale are kept independent. The two-loop 0{at) correction to the Higgs mass from the matching 
between the MSSM and the SM in the EFT approach can be written as the sum of various 
contributions: 

( 12 ) 
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The meaning of the various terms in this equation is explained below. The term is the 

contribution from the effective potential in the DR scheme, which was calculated by Espinosa and 
Zhang [15]: 
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where Xt = At — jj, cot (3, Yt = At + jj, cot /I, C /3 = cos (3, sp = sin /I, we use the notation i = zjmi 
where stands for any of the parameters p, Xt or F^, and the dehnitions 
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Here Lg(x) is the dilogarithm function. Below the SUSY scale we use the SM as an effective held 
theory in the MS scheme. Then we need to write the MSSM top mass and the EW vev (in the 
DR scheme) in terms of the SM ones in the MS scheme in the one-loop 0{at) correction to the 
Higgs mass. Doing so will produce an additional (shift) contribution at two loops 
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Unlike the two-loop 0{atas) correction to the Higgs mass, the 0{a^) one receives a wave-function 
renormalization contribution. It arises as a combination of the one-loop 0{at) contribution of the 
stops to the wave-function renormalization of the Higgs held and the one-loop correction to the 
Higgs mass from the matching at the SUSY scale. It reads: 


XU ln^+X!--Xt 
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Finally, we need to subtract the 0{a^) corrections to the Higgs mass associated with the contribu¬ 
tion of the top-quark loops to the effective potential because it is already present in the matching 
at the EW scale. The two-loop 0{af) correction to the Higgs mass in the SM from the matching 
at the top mass, which receives EP and WFR contributions, is given in eq. (20) of [24]. We extract 
the EP piece which is given by 
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Evaluating eq. (12) we obtain for the Higgs quartic coupling 
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After taking into account all the contributions in eq. (12), we checked that the logarithmic de¬ 
pendence on the top mass of the SM quartic coupling is canceled, as it is shown in eq. (21). We 
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also verified analytically that the inclusion of the two-loop 0{al) correction in eq. (21) makes the 
result of the pole Higgs mass independent of the renormalization scale at this order. 

The two-loop 0{atas) correction was also re-computed in this work. The explicit expressions 
for the SUSY thresholds are too long to be reported here and can be found in the SusyHD 
package [27]. 


On-shell scheme 


A change in the renormalization of the parameters entering in the one-loop SUSY thresholds 
to A will produce a two-loop (shift) contribution. We present the relation between the MSSM 
parameters in the DR and OS schemes. In particular, we need the relations for the stop masses 
and mixing at 0{as) and 0{at), the latter for degenerate stops. This will determine the shift 
contributions to the two-loop 0{atas) and 0{af) SUSY corrections in the OS scheme. 

In the OS renormalization scheme the masses are dehned as the poles of the propagators. The 
relation between the DR and OS masses for a scalar particle with squared mass is given by 

^2 (OS) _ ^2(DR)(g^ _ ^22) 

6m\Q) =Refl{m\Q), (23) 


where (Q) is the tree-level DR mass evaluated at the renormalization scale Q, and n(m^, Q) 

is the DR renormalized one-loop self-energy. 


On the other hand, the OS renormalization for the mixing angle is more subtle. At tree-level, 
the mixing angle of the stops is 
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We use the symmmetric renormalization for the stop mixing angle (for a discussion on possible 
renormalizations see [91] and references therein): 
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where IIi 2 (p^) is the off-diagonal self-energy of the stops. We dehne the OS combination [rritXt)^^ 
from eqs. (24) and (25), which implies 
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In the usual dehnition for the stop mixing on-shell, = (mtXt)^^/, terms proportional to 
\og{mi/mt) appear in the two loop thresholds. In the EFT approach these logs are big and need 
resummation. Therefore we use a different dehnition for Xt which does not produce such terms 
and is more suitable for the EFT computation: 
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where the numerator is computed from eq. (26) and we stress again that is the top mass in 

the SM as any other MS quantities in this paper. An analogous dehnition applies for the sbottom 
and stau mixings. The decoupling of heavy particles like the gluino is ensured in our on-shell 
renormalization scheme. 


For the squarks, the 0{as) shift (neglecting the quark masses) reads [43]: 
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and for the shift between the DR top mass in the MSSM and the MS top mass in the SM (which 
is also given in [26]) 
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The functions Fio and Fu in eq. (29) are dehned as: 
1 — a;^ 


Xio(x) = 


x^ 


Fii{xi, X 2 ) — —2 -|- 


In |1 — a;"^| 

2 (a;^ Ina;^ — x^ \n.x\) a;f(l — x\Y Yi\l — X 2 \ — xl{l — xfY In |1 — a;^ 


^2 _ 

vT 1 vT o 


+ 


_ ™ 2 ™ 2/'™2 _ ^24 

n -^2 •ria'2t-^l X 2 ) 

Analogously for the 0{a‘^) corrections for degenerate scalars [15] 
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As it was discussed in section 3.2, in the DR scheme there are power-like corrections from the 
gluino-stop loops to the Higgs quartic coupling which do not decouple in the limit of heavy gluino. 
We illustrate this effect for the simplihed case of degenerate stops 
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While our on-shell result, obtained from the DR one by shifting the parameters in the one-loop 
0{at) correction, guarantees the decoupling of heavy gluino 
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We also see that eq. (37) does not contain large logarithms Inm^/m*. At last, the two-loop 0{al 
SUSY threshold in our on-shell scheme is given by 
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